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Abstract 

The effects of initial conditions and system parameters on entanglement dynamics and asymptotic 
entanglement for a two-qubit anisotropic XY Heisenberg system in the presence of an inhomo- 



geneous magnetic field and spin-orbit interaction are investigated. We suppose that each qubit 
interacts with a separate thermal reservoir which is held in its own temperature. The effects of the 
parameters of the system and environment like spin-orbit interaction and temperature difference 
'p ■ of reservoirs are also discussed. 
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I. INTRODUCTION 



Entanglement is a central theme in quantum information processing which is first noted 
by Schrodinger It strongly affects our conceptual implication on physics, and force 

us to change significantly our perspective of Nature. Entanglement implies that the best 
knowledge of the whole of a composite system may not include complete knowledge of its 
parts. In mathematical sense a pure state of pair of quantum systems is called entangled 
if it is unfactorizable. A mixed state p of a bipartite system is said to be separable or 
classically correlated if it can be expressed as a convex combination of uncorrelated states 
Pa and ps of each subsystems i.e. p = J2i ^iP % a ® Pb sucn that > and J2i u i — 1j 
otherwise p is entangled js, 4]. Entanglement has no classical analog and then can be con- 
sidered as a uniquely quantum mechanical (non-classical) resource that plays a key role 
in many of the most interesting applications of quantum computation and quantum infor- 
mation processing such as: quantum teleportation, entanglement teleportation, quantum 
cryptography, and etc. 0, 4]. Performance the above mentioned tasks needs to quantifying 
and optimizing the amounts of the entanglements in a suitable multipartite quantum sys- 
tem. Many measures of entanglement have been introduced and analyzed ^, |], 6], but the 
one most relevant to this work is entanglement of formation, which is intended to quantify 
the resources need to create a given entangled state [51]. For the case of a two-qubit sys- 
tem Wootters [5] has shown that entanglement of formation can be obtained explicitly as: 
E(p) = E[C(p)} = h ( 1 +^+ c2 j ( where h(x) = — xlog 2 x — (1 — x) log 2 (l — x) is the binary 
entropy function and C(p) = max{0, 2X max — J2i=i ^i} is the concurrence of the state, where 
AjS are positive square roots of the eigenvalues of the non-Hermitian matrix R = pp, and 
p is defined by p := (cr y (g) a y )p*(a y ® a y ). The function 5 is a monotonically increasing 
function and ranges from to 1 as C goes from to 1, so that one can take the concurrence 
as a measure of entanglement in its own right. In the case that the state of the system is 
pure i.e. p = \ip)(ip\, \ip) = a\00) + b\01) + c|10) + the above formula is simplified to 

C(\ip)) = 2\ad-bc\. 

Real Quantum systems are not isolated from their environment. Unavoidable interaction 
between system and its environment cause to leakage the coherence of the system to envi- 
ronment and hence quantum-to-classical transition (see [7| and references therein) . More 
precisely, this interaction will, in general, create some entanglement between the states of 
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quantum system and the huge states of environment and as a consequence, quantum co- 
herence initially localized within the system will become a shared property of composite 
system-environment state and can no longer be observed at the level of the system, lead- 
ing to decoherence. Decoherence destroys the quantumness of the system and hence will 
decrease the useful entanglement between the parts of the system. A given dynamics for 
a composite quantum system can exhibit several distinct properties for asymptotic entan- 
glement behavior, like: asymptotic death of entanglement (ADE) i.e. entanglement vanishes 
exponentially in time, entanglement sudden death (ESD) i.e. entanglement vanishes faster 
than local coherence of the system, asymptotically steady state entanglement (ASE) i.e the 
entanglement of the system reaches a stationary value, asymptotically, entanglement sudden 
birth (ESB) i.e. an initially separable state acquire entanglement asymptotically and etc 

im. 

Among the numerous concepts considered to implement quantum bits (qubits), approaches 
based on semiconductor quantum dots (QDs) offer the great advantage that ultimately a 
miniaturized version of a quantum computer is feasible. Indeed, Loss and DiViencenzo 
initially proposed a quantum computer protocol based on electron spins trapped in semi- 
conductor QDs in 1998 jiol. [ill. 12]. Here qubit is represented by the spin of a single electron 
in a QD, which can be initialized, manipulated, and read out by extremely sensitive devices. 
Such systems are more scalable and have a longer coherence time than other systems such 
as quantum optical and NMR systems. In this paper quantum entanglement dynamics of 
an open two-qubit system is realized by considering two electrons confined in two coupled 
quantum dots (CQDs) interacting with two separate reservoirs. We refer to this two-qubit 
system as "nanosystem" in the rest of the paper. Because of weak lateral confinement elec- 
trons can tunnel from one dot to the other and spin-spin and spin-orbit interactions between 
the two qubits exist. Thus we can model this nanosystem by a two-qubit spin chain includ- 
ing spin-orbit interaction. The spin chain is the natural candidates for the realization of 
entanglement and Heisenberg model is the simplest method for studying and investigating 
the behavior of the spin chains. On the other hand, in what follows we model the envi- 
ronment by a thermal reservoir and assume that the inter dot separation is large enough 
such that each dot couples to a separate thermal reservoir (bosonic bath). Here the bathes 
are assumed to be in thermodynamical equilibrium at different temperature = ^jr- In 
general, there are two different ways for connecting the quantum dots to the reservoirs: (i) 
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" direct geometry" ; where a high temperature bath is in contact with a QD in the presence 
of a strong magnetic field i.e. bAT > and (ii) "indirect geometry" ; where a high temper- 
ature bath is in contact with a QD in the presence of a weak magnetic field i.e. bAT < 0. 
The results show that the non-equilibrium thermal entanglement dynamics depends on the 
geometry of the system. 

Entanglement properties of Heisenberg systems at thermal equilibrium (thermal entangle- 
ment) are extensively studied after Nielsen [ljj, who first studied the entanglement of a 
two-qubit Heisenberg XXX chain (see [lij] and references therein) . For non-equilibrium 
thermal entanglement in spin systems, Eisler et al. [15| calculated the von Neumann en- 
tropy of a block of spins in XX spin chain in the presence of the energy current and showed 
that the enhancement of the amount of entanglement due to an energy current is possible. 
After them, the non-equilibrium therma 
been studied in a number of works 17 



entanglement for steady state of some systems has 

nrx n 

181 . |19| ]. For example, Quiroga in Ref. [20| consid- 



ered a simple spin chain system (XXX-Heisenberg) which is in contact with two different 
heat reservoirs and showed that for the steady state, a temperature gradient can increase 
or decrease entanglement depending on the internal coupling strength between spins. Dy- 
namics of non-equilibrium thermal entanglement of the same system has been studied by 
Sinaysky et al. 21[. They have derived an analytical expression for the density matrix of 
the system at a finite time. They also have shown that the system converges to a steady 
state, asymptotically and the amount of entanglement of the steady state takes its maxi- 
mum value for unequal bath temperatures and also the local energy levels can maintain the 
entanglement at higher temperatures. However, the non-equilibrium thermal entanglement 
dynamics of more involved spin systems (e,g XY and XYZ Heisenberg systems) has not 
been considered, yet. In this paper, we will investigate dynamics of non-equilibrium ther- 
mal entanglement of a two qubit anisotropic XY Heisenberg system (nanosystem) in the 
presence of the inhomogeneous magnetic field and the spin orbit interaction. The influence 
of the parameters of the nanosystem (i.e. magnetic field (B), inhomogeneity of magnetic 
field(b), partial anisotropy(x), mean coupling (J) and the spin-orbit interaction parameter 
(D)) and environmental parameters (i.e. temperatures 7\ and T 2 or equally Tm and AT, and 
the couplings 71 and 72) on the entanglement of the nanosystem is investigated. We have 
shown that, there is a steady state entanglement for asymptotically large times. The size of 
this steady state (asymptotic) entanglement and the dynamical behavior of entanglement 
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depend on the parameters of the model and also on the geometry of the system. Increasing 
the temperature difference AT, and mean temperature Tm, decrease the amount of asymp- 
totic entanglement. We have also shown that, the size of T£p (temperature over which the 
entanglement vanishes) and the amount of entanglement can be improved by adjusting the 
value of the spin-orbit interaction parameter D. The maximum entanglement (C = 1) can 
be achieved for the case of large values of D and zero temperature reservoirs(Ti = T2 = 0). 
Furthermore, we find that the indirect geometry is more suitable for creating and main- 
taining the entanglement. The results obtained here are consistent with those obtained in 



14. 20. 



2l|. 



The paper is organized as follows. In Sec. II we introduce the Hamiltonian of the whole 
system-reservoir under the rotating wave approximation and then write the Markovian mas- 
ter equation governed on the nanosystem by tracing out the reservoirs' degrees of freedom. 
Ultimately, given some initial states, the density matrix of the nanosystem at a later time 
is derived exactly. The effects of initial conditions and system parameters on the dynamics 
of entanglement and entanglement of asymptotic state of the nanosystem are presented in 
Sec. III. Finally in Sec. IV a discussion concludes the paper. 

II. THE MODEL AND THE HAMILTONIAN 

The total Hamiltonian of the nanosystem which is interacting with two heat reservoirs 
(bosonic bath) is described by 

H = Hg + Hbi + H B 2 + H S B1 + HsB2, (1) 

where H$ is the Hamiltonian of the nanosystem, Hbj is the Hamiltonian of the jth bath 
(j=l,2) and HsBj denotes the system-bath interaction Hamiltonian. Nanosystem consists of 
two spin electron confined in a two coupled quantum dots and is described by a two-qubit 
anisotropic Heisenber g X Y-model in the presence of an inhomogeneous magnetic field and 



spin-orbit interaction 14j with the following Hamiltonian 

H s = \{J X *{o* + J y a\a\ + B x ■ a x + B 2 a 2 

+ D • (<r x x a 2 ) +8 a x Y ■ cr 2 ) (2) 

where o-j = (crj, crj, cj) is the vector of Pauli matrices, Bj (J = 1, 2) is the magnetic field on 
site j, (/i = x, y) are the real coupling coefficients (the chain is anti-ferromagnetic (AFM) 



for J M > and ferromagnetic (FM) for < 0) and D is Dzyaloshinski-Moriya vector, which 



is of first order in spin-orbit coupling and is proportional to the coupling coefficients 
r is a symmetric tensor which is of second order in spin-orbit coupling 



22 



23 



24 



) and 



25|]. For 



simplicity we assume Bj = Bj z such that B\ = B + b and B2 = B — b, where b indicates 
the amount of inhomogeneity of magnetic field. The vector D and the parameter 5 are 
dimensionless, in system like coupled GaAs quantum dots |D| is of order of a few percent, 
while the order of last term is 10~ 4 and is negligible. If we take D = JD z and ignore the 
second order spin-orbit coupling, then the above Hamiltonian can be written as: 



H s = JxiafaZ + °"i °2 ) + ^(1 + iB)o\o 2 + J(l - iD)a l o\ 
,B + b. . ,B-b. z 
+ (—3- K + (—3— (3) 

where J := Ja: l" Ja , is the mean coupling coefficient in the XY-plane, y := , specifies the 

* <J x\ <Jy 

amount of anisotropy in the XY-plane (partial anisotropy, — 1 < x < 1) and a ± = ^(cr x ±i(T y ) 
are lowering and raising operators. The spectrum of Hs is easily obtained as 

|e 1)2 ) = |*±> = iV±((-^)|01) + |10» , e 1 , 2 = ±e , 

(4) 

\e 3A ) = |£±) = M±((^)|00) + |11)) , 53,4 = ±V- 

Here the eigenstates are expressed in the standard basis {|00), |01), |10), |11)}. In the above 
equations N± = (1 + ji+^p )~ 1/2 and M ± = (1 + (^J 2 ) 2 ) -172 are the normalization 
constants. Here we define, £ := (6 2 + J 2 + (JD) 2 ) 1 / 2 and 77 := (£? 2 + (Jx) 2 ) 1//2 , for later 
convenience. 

The Hamiltonian of the reservoirs for each spin j = 1 , 2 are given by 

H B j = J2 u n b ij b nj- ( 5 ) 

n 

The interaction between the nanosystem and the bosonic bathes in the rotating wave ap- 
proximation is as following 

H SBj = 4 E rf^i + ^7 E = E v <■<■■(,!, (6) 

n n ^ 

The system operator Vj ;fM are chosen to satisfy [Hs,Vj !fi ] = tuj^Vj^, and the fj s are the 
random operators of reservoirs and act on the bath degrees of freedom. Physically, the 



index fi corresponds to transitions between the internal levels of the nanosystem induced 
by the bath. The dynamics of the whole nanosystem+reservoirs is described by a density 
operator (a) satisfying the Liouville equation a = — i[H, a]. If the coupling strengths of 
the nanosystem and the environment are weak, the evolution of the nanosystem does not 
influence the states of the reservoirs and one can write a(t) = p{t)pBi (0)pb2(0) (irreversibility 
hypothesis), where p(t) is the reduced density matrix describing the nanosystem and each 
bosonic bath is described by a canonical density matrix of the form pBj = e~^ jBBi jZ , where 
Z = Tr(e _/ ^ H£y ) is the partition function of the jth bath. 

In the Born-Markov approximation the master equation describing the dynamics of the 



reduced density matrix of the nanosystem is [16|, |26] 

dp 



(jf - i [H s ,p] + £ 1 (p) + £ 2 (p) (7) 



where £j(p) are dissipators given by [26j 



£Ap) = E j £K^>){K,m> \vLpW - (i - ^ s '~)\yj,»v},M- (8) 

fl, V 

here J^l(^j,u) is the spectral density of the jth reservoir, 




where G al 3(r) is the environment self-correlation function, 

Gap(r) =Tr Bj [p Bj f ju (T)f jai ] (10) 

and fj u (t) = e~ lHB i T f^ fl e lHB: > T . Spectral densities encapsulate the physical properties of 
the environment and play an immensely important role in the theoretical and experimental 
studies of the decoherence. In this paper, we will consider the bosonic thermal bath as 
an infinite set of harmonic oscillators and apply a Weisskpof-Wignner-like expression for 
spectral density such as J^{oj^) = ^jioj^n^uj^ where rij{u^) = (eP^ — 1) _1 denotes the 
thermal mean value of the number of excitations in the jth reservoir at frequency uj^ and 
temperature Tj = j- and 7j(^ M ) is the coupling strength of nanosystem and the jth reservoir. 
For simplicity we take 7j(^ M ) = jj- Thus, the dissipators £j(p) become 
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+ E jW W( 2 ^,J^i,/*-{p.^,AJ} + )) (ii) 



with the transition frequencies 



and the transition operators 



^2 + uj 3 = -uj 2 , (12) 

V},2 = a i)2 |* + )(E-|, (13) 
T>- 3 = a,, 3 |^)(S + |, 
t>, 4 = aj-,4|*">(E"|, 

(14) 



where 



a h 2 | 2 = I a Jl3 | 2 = — (^ - J 2 X - (~l) j Bb), (15) 



1 

can be obtained from the spectrum of the nanosystem Hamiltonian. 



The quantum master equation (j7]) has an important property, when the spectrum of H s 
(see eq.(plj) is non-degenerate: In the energy basis, {\Ei)}f =1 , the equations for diagonal 



elements decouple from nondiagonal ones 16j. Furthermore, nondiagonal elements are not 



coupled and the time dependence of these elements has the simple form 

(>;.,(/) !>,,,(())("■'. (16) 

where Ojj G C are determined by the nanosystem parameters. The equations for diagonal 
elements have the following form 

R(t) = BR(t), (17) 



where dot denotes the time derivative, R(t) = (pn(t) , p 22 {t) , p 33 {t) , p44(t)) T and B is the 
time independent 4x4 matrix 



B 



where 





Y 2 - 



- W + 
Y 2 + 



xt 

Y 2 

{Xt+Y 2 




Y 2 + 

XT 
o 



^ = 2E J®te»A i I 2 , 



5=1,2 



The analytical solution of the equation (fTTI) in the energy basis is given by 



(18) 



(19) 



R(t) = M(t)R(0) 



(20) 



where M(t) = [ ?7ijj]4x4, and the elements uiij are defined by 
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mi i = {Xt + Xr e"^ ) (y + + Y 2 ] 

-^1*2 

m l2 = -^ V (l-e- tx ^(l-e- tY *)X+Y+, 
1 'l-e-^)X+(Y+ + Y 2 -e-^), 



m 42 



m 14 = (X+ + Xre-' Xl )(l-e-' y2 )r 2 -, 



X 4 r 2 
1 



™2i = -^rril ~ e- tx ^){l - e-^)X-Y 2 , 

^1*2 

m 2 2 = -\ r (xr + x 1 + e-^)(r 2 - + r 2 + e -^; 

™23 = Y V^ +X l +e ^ 1 )( 1 - e_<y2 ^ ; 
-M-T2 

^24 = ^V(i-e- tXl )xr(y 2 - + y+ e -^), 

msi = ^V(l-e- tXl )Xr(y 2 + + F 2 -e-^), 

™32 = Y V^ +X l +e ^ 1 ^ 1 - e_<y2 ^ + - 

™ 3 3 = ^ ^ + ^ + 6 ~ tXl ^ + + Y * < 
A1X2 

^34 = -\ r (i-e^)(i-^ y2 )xry 2 + , 
1 

X\~Y 2 

1 l-e-^)x+(y- + r+ e -^), 



m 14 = (X^ + Xfe-^Xl-e-^KT, 



2 



^43 = ^(l-e-^Xl-e-^X+Y,- 

m 44 = _L r (JC+ + Xre-«0(i7 + ^-« r »), 
A1X2 



(21) 



here we have defined: X^ = + X^ and = K+ + . 

There is a singular point £ = r], for which the spectrum (j3J) becomes degenerate and the 
above solution is not valid. The state of the system is not well defined at this critical 
point. This critical point assigns a critical value for the parameters of the nanosystem such 
as critical magnetic field (B c ), critical parameter of inhomogeneity of magnetic field (6 C ), 
critical spin-orbit interaction parameter (D c ) and etc. The behavior of the entanglement 
of the system changes abruptly when the parameters cross their critical values (see section 
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III). 

In the following we will examine a class of bipartite density matrices having the following 
standard form as the initial state of the system {27] 



Ps(0) 



(22) 



p+ v 

W! z 

z* w 2 

. V yU_ j 

These kind of density matrix are called X states class and arises naturally in a wide variety 
of physical situations. This class contains some important subsets like the pure Bell states, 
the states which can be expressed as a mixture of Bell states, Werner states and so on. If 
the initial state, p s (0), belongs to the set of X states f[2"2"l) then Eq (TJJ guarantees that p s (t) 
given by Eqs. (I2"3l and [To]) also belongs to the same set. Therefore, the only non-vanishing 
off diagonal components of the density matrix in the energy basis are 



Pn(t) = p 12 (0)e- m - t{Xl+Y2)/2 , 



P2i(t) = Pl2 (ty, 

P43<» = P34(t)*. 



(23) 



Knowing the density matrix we can calculate the concurrence C(p(t)) = max{0, 2\ max (t) 
E- =1 \(t)} where, 



Al, 2 (*) = I \fp S n(t)Ps44(t)± I Psu(t) 



As,4(*) = I \f Ps22(t)ps33(t)± I ps23(t) 



(24) 



Unfortunately, the Aj(t)s depend on the parameters involved. This prevents us from 
writing an analytical expression for concurrence. But it is possible to evaluate concurrence, 
numerically for a given set of the parameters. The results are shown in Figs. 1-7. Figs. 
1-3 depict the dynamical behavior of concurrence versus parameters of the nanosystem 
and reservoirs and Figs. 4-7 illustrate the asymptotical behavior of the concurrence versus 
parameters of the nanosystem and reservoirs. Without loss of generality we can assume 
that J > 0, since the above formula are invariant under substitution J — > — J. This means 
that the dynamical behavior of the FM chain is the same as the AFM chain. 
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A. Asymptotic case 



For a class of states, p s t, the dissipative and decoherence mechanisms (second and third 
terms in the master equation (j7j)) compensate the unitary dynamics which is governed by 
nanosystem Hamiltonian (first term in master equation (j7j)) i.e. i[Hs, Pst] = £i(pst) + £2(f>st) 
or jj^pst = 0. These states are called stationary states because they are constant in time. 
If there exist such a stationary state solution for the master equation, the system tends 
to it asymptotically in large time limit i.e. lim^oo p(t) — > p aS ym. = Pst- For the present 
system (thermal reservoirs and interacting nanosystem), nondiagonal elements ffT6j) vanish 
asymptotically at large time limit and hence p(t) converges to a diagonal density matrix (in 
the energy basis) with elements not depending on the initial conditions: 

1 



Pasym. 



x x y 2 



diagonal(X+Y+, X^Yf, X^Y 2 + , X+Yf). 



(25) 



The asymptotic concurrence is given by C(p, 



asym. 



C°° = max{0, 2\ max - £■=! A,} with 



asym. asym. • i asym. 
Psll PsAA ± I PslA 



A 



3, 1 



asym. asym. i i asym. 
Ps22 Ps33 ^ I Ps23 



(26) 



where 



asymp. 
Psll 



asymp. 
Psl4 



asymp. 
Ps22 



asymp. 
Ps23 



asymp. 
Ps33 



asymp. 
Ps44 



1 



2r t X 1 Y 2 

Jx 

2r ] X 1 Y 2 
1 

2£XY 2 
J(l + iD 



(( V + B)X^Y 2 + + ( V -B)X+Y 



2 h 



(x 1 -F+-x 1 + r 2 -)= P -r p - 



((^ + 6)x 1 + y+ + (^-6)x 1 -y 2 



1 



■(x ] f y 2 + — x l y 2 



(p. 



2 h 



asymp. \ * 



*32 



2^y 2 
i 



((e-&)x+y 2 + + (e + 6)xry 2 



((7 7 - J B)xry 2 + + (7 7 + J B)x 1 + yr 



2 V) 



(27) 



2/7X^2 

There is an interesting limiting case for which the coupled QDs are in contact with two 
independent reservoirs at identical temperatures (Pi = p 2 = P). In this case, it is easy to 
show that 
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Y + 
Y 2 



e ^wi _|_ i ■ 



*2 



e 0wi + 1 

1 



e/awa + i ' y 2 e /^ 2 + i 

By substituting these relations into the Eq. fl25l) . the reduced density matrix p asym , takes 
the thermodynamical canonical form for a system described by the Hamiltonian H$ at 
temperature T = as expected. This means that 

e -PHs 



pasymX^T = 0) = p T 



(28) 



where Z = Tr(e l3Hs ) is the partition function. Thermal entanglement properties of such 
systems have bee,, studied substantially, iu out previous work Q. Thus, for the special 
case AT = 0, the results coincide with the results of Ref.|14|. 



III. RESULTS 



The non-equilibrium thermal concurrence as a function of time for three values of tem- 
perature difference (AT = T\ — T 2 ) and for a fixed value of mean temperature (Tm = Tl ^ T2 ) 
are plotted in Fig. 1, for the case of "direct geometry" of connection (6AT > 0). The pres- 
ence of temperature difference has not effective influence on the dynamics of entanglement 
at early times of evolution, but changes effectively the behavior of the asymptotic entangle- 
ment (which is more evident form Figs. 4-7). Fig. 2 depicts the variation of entanglement 
dynamics for some values of Tm and for a fixed AT in the case of "direct geometry" of 
connection (bAT > 0). By Increasing Tm, thermal fluctuations suppress quantum fluctua- 
tions and hence decrease coherent oscillations at early times of evolution and also decrease 
asymptotic entanglement. 

Perhaps surprisingly, the decoherence due to environmental interaction does not prevent the 
creation of a steady state level of entanglement, regardless of the initial state of the system. 
This is demonstrated in figure 3 which shows the time evolution of the non-equilibrium 
thermal concurrence for a given set of parameters and for four different initial states: (i) a 
maximally entangled state, the Bell state, |^(0)) = ^(|01) + 1 10) ) (ii) a separable state, 
| V (0) ) = 1 01) (iii) a mixed state, defined as an equal mixture of a Bell state and a product 
state, e.g. p s (0) = |(|00> + |11>)(<00| + (11|) + ||01><01| and finally (iv) an unpolarized 
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state, p s (0) = \l. Despite the presence of decoherence (due to interaction with environ- 
ment) the results of figure 3 show that the concurrence reaches the same steady state value, 
C°° (after some oscillatory behavior) for a given set of parameters, regardless of initial state 
of the system. Clearly the Heisenberg interaction in the Eq. ([3]) serves to maintain an 
entangled asymptotic state despite the presence of decoherence. At early times of evolution 
the amount of D determines the frequency of oscillations. Each plot in Fig. 3 contains two 
cases a) D < D c , in this case asymptotic value of entanglement decreases with D and b) 
D > D c , in this case asymptotic value of entanglement increases as D increases. 
In Fig. 4, the asymptotic non-equilibrium thermal concurrence is plotted versus Tm and D 
for different values of AT. For the case of identical temperatures (AT = 0) the results are 
the same as ref. [141 ] . as mentioned in the previous section. This figure shows that there is a 
critical mean temperature (T^f) over which entanglement vanishes (ESD phenomenon). The 
size of Tm and the amount of entanglement can be improved by increasing D. For Small 
values of D (D < D c ), an increase in AT increases the size of T^ i.e entanglement can 
exists in higher mean temperatures due to existence of temperature difference of reservoirs 
(environment induced entanglement). 

The variation of the asymptotic non-equilibrium thermal concurrence as a function of AT 
and D for fixed values of Tm and b is illustrated in Fig. 5. The behavior of concurrence 
is dependent on the geometry of connection. For the case of "direct geometry" (bAT > 0) 
and for D < D c , no entanglement is observed but for the same geometry and for D > D c 
the amount of entanglement is nonzero and it increases as D or AT increases. For the case 
of "indirect geometry" (6AT < 0), there is a nonzero entanglement for all values of D. The 
amount of entanglement is an increasing function of D and increases with AT for the values 
of AT < and decreases with AT for the values of < AT < Tm- 
Figs. 6 and 7 show the behavior of the asymptotic non-equilibrium thermal concurrence ver- 
sus Tm and AT for different values of D and for the symmetric (b = 0) and nonsymmetric 
{b 7^ 0) cases, respectively. Both figures reveal that, increasing D cause to the appearance 
of entanglement in the larger region of Tm — AT plane. The departure between symmetric 
(Fig. 6) and nonsymmetric (Fig. 7) case is more obvious for D < D c . For nonsymmetric 
case, Fig. 7 shows that the "indirect geometry" of connection is more suitable for creating 
entanglement when D < D c . In both symmetric and nonsymmetric cases, maximum entan- 
glement (C = 1) can be achieved in the case of identical temperatures (AT = 0), zero mean 
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temperature (Tm = 0) (i.e when both of reservoirs are in the ground state (Tj. = T<i = 0)) 
and large values of D (it is in agreement with the results of 14J]). 



IV. DISCUSSION 



The Dynamics of non-equilibrium thermal entanglement of an open two-qubit nanosys- 
tem is investigated. The inter-qubit interaction is considered as the Heisenberg interaction 
in the presence of inhomogeneous magnetic field and spin-orbit interaction (arises from the 
Dzyaloshinski- Moriya (DM) anisotropic antisymmetric interaction). Each qubit interacts 
with separate thermal reservoir (bosonic bath) which is held in its own temperature. The 
effects of the parameters of the model, including the parameters of the nanosystem (espe- 
cially, the parameter of the spin orbit interaction (D)) and environmental parameters (par- 
ticularly, mean temperature Tm and temperature difference AT)), on the non-equilibrium 
thermal entanglement dynamics of the nanosystem is investigated, by solving the quantum 
Markov-Born master equation of the nanosystem. An analytical solution of the master 
equation is derived and then entanglement dynamics and asymptotic entanglement of the 
nanosystem versus the parameters of the model is studied. Resolving the entanglement 
dynamics allowed us to distinguish between entanglement induced by the interaction and 
by the environment. The results show that, decoherence induced by thermal bathes are 
competing with inter-qubit interaction terms to create a steady state level of entanglement, 
as measured by the concurrence. The size of this steady state (asymptotic) entanglement 
and the dynamical behavior of entanglement are dependent on the parameters of the model 
and also depend on the geometry of connection. Increasing temperature difference AT, 
and mean temperature Tm, kill the asymptotic entanglement. Indeed, thermal fluctuations 
suppress the quantum fluctuations (i.e. all quantum effects such as the entanglement and 
local coherence), and hence the entanglement of the system dies at a critical temperature 
T^' (ESD). We have shown that, the size of T^ and the amount of entanglement can be 
enhanced by choosing a suitable value of spin-orbit interaction parameter D. The maximum 
entanglement (C = 1) can be achieved for the case of large values of D and zero temperature 
reservoirs(Ti = T2 = 0). For physical realization of the model we address a two coupled 
quantum dots which are interacting with two independent thermal bathes. Furthermore, we 
find that the indirect geometry of connection is more suitable for creating and maintaining 
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the entanglement. In this case, the heat current between two QDs is substantially decreased 
but the concurrence can enhances by temperature difference AT | 20l| . For the special case 
AT = 0, our results confirm the results of 14j. The results of [20|, [21( are also obtained 
for the special case of D = \ = and by considering b = e\ — e 2 . The results can provide 
useful recipes for realistic quantum information processing in noisy and non-equilibrium 
environments. 
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FIG. 1: Dynamics of non-equilibrium concurrence for the initial reduced density matrix p s (0) = 

^(|01)(01| + |10)(10|). The parameters of the model are chosen to be 7l = 72 = 0.02, J = 1, 

X = 0.9, B = 2, b = 1, Tm = 1.5 and for different values of temperature difference AT: (a) AT = 

(b) AT = 1 (c) AT = 2. All parameters are dimensionless. 
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FIG. 2: Dynamics of non-equilibrium concurrence for the initial reduced density matrix p s (0) = 

^(|01)(01| + |10)(10|). The parameters of the model are chosen to be 7l = 72 = 0.02, J = 1, 

X = 0.9, B = 2, b = 1, AT = 0.5 and for different values of mean temperature Tm- (a) Tm = 1 

(b) Tm = 1.5 (c) Tm = 2. All parameters are dimensionless. 
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FIG. 3: (Color online) Dynamics of non-equilibrium concurrence for different values of D and 
different initial state. The parameters of the model are chosen to be 71 = 72 = 0.02, J = 1, 
X = 0.9, B = 2, b = 0.5, Tm = 1 and AT = 0.5. Each plot contains two graphs for (a) D < D c (b) 
D > D c (D c ~ 1.8868). All parameters are dimensionless. 
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FIG. 4: (Color online) Asymptotic entanglement vs. Tm and D, The parameters of the model are 
chosen to be 7l = 72 = 0.02, B = 4, J = 1, x = 0-3 and b = -3.5 for (a) AT = (b) AT = 1 (c) 
AT = 2 (d) AT = 3 {D c ~ 1.68523). All parameters are dimensionless. 




FIG. 5: (Color online) Asymptotic entanglement vs. AT and D. The parameters of the model are 
chosen to be 7l = 72 = 0.02, B = 4, J = 1, X = 0.3 and b = -3.5 and T M = 2 (D c ~ 1.68523). All 
parameters are dimensionless. 
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FIG. 6: (Color online) Asymptotic entanglement vs. Tm and AT. The parameters of the 
are chosen to be 7l = 72 = 0.02, B = 4, J = 1, x = 0-3 and b = for (a) D = (b) TJ = 
D f» Z) c (d) D = 5 {D c ~ 3.88458). All parameters are dimensionless. 
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FIG. 7: (Color online) Asymptotic entanglement vs. Tm and AT. The parameters of the 
are chosen to be 7l = 72 = 0.02, B = 4, J = 1, x = 0-3 and b = 1 for (a) D = (b) TJ = 
-D f» Z) c (d) D = 5 (-D c ~ 3.75366). All parameters are dimensionless. 
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